Critical dynamics at incommensurate phase transitions and NMR relaxation 

experiments 
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We study the critical dynamics of crystals which undergo a second-order phase transition from a 
high-temperature normal phase to a structurally incommensurate (/C) modulated phase. We give 
a comprehensive description of the critical dynamics of such systems, e.g. valid for crystals of the 
A2BX4 family. Using an extended renormalization scheme, we present a framework in which we 
analyze the phases above and below the critical temperature Tj. Above Tj, the crossover from the 
critical behavior to the mean-field regime is studied. Specifically, the resulting width of the critical 
region is investigated. In the IC modulated phase, we consider explicitly the coupling of the order 
parameter modes to one- loop order. Here the Goldstone anomalies and their effect on measurable 
quantities are investigated. We show their relation with the postulated phason gap. While the theory 
can be applied to a variety of experiments, we concentrate on quadrupole-perturbed nuclear magnetic 
resonance (NMR) experiments. We find excellent agreement with these dynamical measurements 
and provide answers for some questions that arose from recent results. 



I. INTRODUCTION 



This paper is concerned with the critical dynamics 
of crystals undergoing a second-order phase transition 
from a high-temperature normal (iV) phase to a struc- 
turally incommensurate (IC) modulated phase. In the 
IC phase, the translational symmetry of the lattice is 
broken by a modulation in such a way that the char- 
acteristic wave vector q/ is a non-rational multiple of a 
basic lattice vector. The occurrence of incommensurate 
modulations is in general-understood as a consequence 
of competing interactions .□ The most important charac- 
teristic of these systems is that the ground state does 
not depend on the actual phase of the incommensurate 
modulation at a lattice site. This implies that the initial 
phase of the modulation wave is arbitrary and one must 
take into account a phase shift degeneration of the ground 
state energy. Consequently, not only the amplitude of the 
modulating vector is required to characterize each config- 
uration, but in addition the phase at a arbitrary lattice 
site must be fixed. Therefore, a two-dimensional order 
parameter (OP) has to be employed in order to describe 
the phase transition from a N phase to an IC modulated 
phaseH Interesting static properties, e.g., the very rich 
phase diagrams of systems with competing interactions, 
emerge B 

However, in this work we concentrate on dynamical 
properties. Considering fluctuations of the OP, the 
normal modes can be expressed in terms of the trans- 
verse and longitudinal coroponents V''" ^iid t/"" in the 
two-dimensional OP space.cl The fluctuations of -0^ and 
can be identified with the fluctuations of the phascj 
and the amplitude of the modulation in the crystal.u 
As a consequence of the OP being two-dimensional, the 



lattice dynamics of structurally incommensurate phases 
shows some peculiar effects which are different from or- 
dinary crystalline phases. Namely, below the transition 
temperature T/ two non-degenerate, branches of modes 
appear in the dynamical spectrum. □ The "amplitudon" 
branch, connected with the fluctuations of the amplitude 
of the incommensurate modulation, exhibits common 
soft-mode behavior. In addition the "phason" branch 
represents the massless Goldstone modes of the system, 
here originating of the invariance of the crystal energy 
with respect to a,ghase shift. Because of the massless 
Goldstone modesoO present in the entire IC phase, new 
types of anomalies may occur. Examples oLsueh anoma- 
lies were discussed in the literature before.Lrii3 Thus we 
expect some peculiar features of the dynamics in the IC 
modulated phase stemming from the Goldstone modes 
and their coupling to the other OP modes. 

The purpose of this paper on the one hand is to pro- 
vide a general framework for the analysis of the criti- 
cal dynamics above and below the N /IC phase tran- 
sition. The theoretical description of such systems is 
based on an 0{2) symmetric time-dependent Ginzburg- 
Landau model, with purely relapcational behavior of the 
non-conserved order parameterB The more general 0{n)- 
symmetric model has been widely studied above the crit- 
ical temperature by means of the dynamical renormal- 
ization group.liirEj Below the critical temperature, the 
0{n) symmetry is spontaneously broken; as mentioned 
above parallel and perpendicular fluctuations have to be 
distinguished. We will start from the fleld-theoretical 
model of incommensurate phase transitions and de- 
rive the corresponding dynamical Janssen-De Dominicis 
functional JIj'Ej which provides us with the framework to 
calculate interesting theoretical properties and correla- 
tions functions, which are required for the interpretation 
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of experimental data. 

Furthermore, we intend to give a comprehensive the- 
oretical description that goes beyond the mean-field or 
quasi-harmonic approach for the IC phase and which is 
missing to date. We present an explicit renormalization- 
group analysis to one-loop order above and below the 
critical temperature T/. The renormalization group the- 
ory will lead us beyond the mean-field picture and pro- 
vide some new insight on the effects of the Goldstone 
modes on the dynamical properties below T/. Some spe- 
cific features of the Gpklstone modes were discussed for 
the statics by LawrieE-3 and for the dynamics in Refs. 
|l| and In the present paper, we extend the analy- 
sis of the 0(n)-symmetric model, specifically for the case 
n — 2. We also consider the crossover behavior from 
the classical critical exponents to the non-classical ones 
in detail, both above and below the critical temperature 
(see also Ref. pj| ). 

Furthermore our model is employed to analyze spe- 
cific experiments. Quadrupole-perturbed nuclear mag- 
netic resonance [NMR) is an established method to in- 
vestigate IC phases in a very accurate way.tS In this 
probe, the interaction of the nuclear quadrupole moment 
(Q) of the nucleus under investigation with the electric 
field gradient (EFG) at its lattice site is measured. The 
fluctuations of the normal modes give rise to a fluctuat- 
ing EFG, which is related to the transition probabilities 
between the nuclear spin levels. As a consequence, the 
relaxation rate 1 /Ti of the spin-lattice relaxation is given 
by the spectral density of the EFG fluctuations at the 
Larmor frequency. We calculate the NMR relaxation 
rate with our theoretical model, and compare our find- 
ings with the experimental data. Our results may be used 
to interpret a variety of experimental findings; however, 
here we will restrict ourselves to the analysis of NMR 
experiments. The theory presented here is appropriate 
for the universality class containing, e.g., the crystals of 
the A2BX4 family. Some very precise NMR experirnenta 
on these crystals were performed over the past years. t3't£l 
Above Tj, these data can be used to analyze the critical 
dynamics in a temperature range ofT — Tj — lOOK more 
closely. Below T/, an identification of relaxation rates, 
caused by fluctuations of the amplitude and the phase, 
respectively, at special points of the NMR frequency dis- 
tribution is possible. Therefore the relaxation rates 1/Ta 
and l/T^, referring to the critical dynamics of the two 
distinct excitations ( "amplitudons" and "phasons"), can 
be studied separately. These experiments led to some 
additional open questions. Above the critical tempera- 
ture T/, a large region was reparted where non-classical 
critical exponents were found.Ea Below T/, the presence 
of a phason gap is discussed in order to clarify some_ex« 
periments as well as the theoretical understandingO'E2l 
We will show how these questions can be resolved within 
the framework of our theory. 

This paper is organized as follows: In the following 
section we introduce the model free energy for a system 
that reveals a to IC phase transition. The dynam- 



ics of the amplitude and phase modes is described by 
Langevin-type stochastic equations of motion, and we 
give a brief outline of the general dynamical perturbation 
theory. In Sec. Ill, the connection between the NMR 
experiments and the susceptibility calculated within our 
theory is discussed. We shall see that the spectral den- 
sity functions are closely related to the measured relax- 
ation times. The high-temperature phase is analyzed in 
Sec. IV. Above Tj, scaling arguments are used to de- 
rive the critical exponents for the relaxation rate. The 
crossover from non-classical to classical critical behavior 
is discussed by means of a renormalization group analy- 
sis, and we comment on the width of the critical region. 
In Sec. V, we apply the renormalization group formal- 
ism to the low-temperature phase. The susceptibility, 
containing the critical dynamical behavior for the ampli- 
tudon and phason modes, is calculated to one-loop order. 
Specifically, the influence of the Goldstone mode is inves- 
tigated. In the final section we shall discuss our results 
and give some conclusions. 

II. MODEL AND DYNAMICAL PERTURBATION 
THEORY 

A. Structurally incommensurate systems 

We want to study second-order phase transitions from 
a high-temperature normal (N) phase to a structurally 
incommensurate (IC) modulated phase at the critical 
temperature Tj. The real-space displacement field corre- 
sponding to the one-dimensional incommensurate mod- 
ulation canJpe represented by its normal mode coordi- 
nates Q(q).cl We will treat systems with a star of soft 
modesLJ consisting only of two wavevectors q/ and — 
along one of the principal directions of tha_ Brillouin 
zone, e.g. substances of the A2BX4 family.El Because 
the incommensurate modulation wave is in most cases, 
at least close to T/, a single harmonic function of space, 
the primary Fourier components (Q(q)) cx 5(q± q/)e*"^° 
with the incommensurate wavevectors ibq/ are domi- 
nating. Using (5(q) as a primary order parameter of 
the normal-to-incommensurate phase transition in the 
Landau-Ginzburg-jWilson free energy functional, diago- 
nalization leads toci 
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X C(ki)Vo (k2)V'f(k3)Vf(k4) SiY^^i) , 

1=1 

with new Fourier coordinates "00 (k) and ijj^Ck) in the 
OP space. Here, we have introduced the abbreviations 
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Below the phase transition, the fluctuations of ipo and 
^|J^ can be identifled with the fluctuations of the phase 
and the amphtude ofJ;he displacement field, named pha- 
son and amplitudonu The wavevector k indicates the 
derivation from the incommensurate wavevector qi. 



(2.2) 



The parameter Tq is proportional to the distance from 
the mean-field critical temperature To/ 



To cx r - To 



(2.3) 



and the positive coupling Uo gives the strength of the 
isotropic anharmonic term. Unrenormalized quantities 
are denoted by the suffix o. 

The functional H[ip"] describes the statics of the 
normal-to-incommensurate phase transition. It repre- 
sents the rt-component isotropic Heisenberg model; in the 
case n = 2, it is also referred to as the XY model. For 
the sake of generality, the n-component order parameter 
case will be considered in the theoretical treatment. 



B. Critical Dynamics 

The critical dynamics of the system under consider- 
ation is characterized by a set of generalized Langevin 
equations for the "slow" variables, which in our case con- 
sist of the non-conserved order parameter fluctuations 
(because of the critical slowing-down in the vicinity oij-a. 
phase transition) .a The purely relaxational behaviorcfl 
is described by the following Langevin-type equation of 
motion 



d_ 
dt 



(2.4) 



The damping is caused by the fast degrees of freedom, 
which are subsumed in the fluctuating forces C". Ac- 
cording to the classification Halperin and Hohenberg, 
we are considering model A£j 

The probability distribution for the stochastic forces 
is assumed to be Gaussian. Therefore 

{C{k,t))=Q, (2.5) 
(r (k, t) C^{k', t')) = 2 Ao (5(k - k') 6{t - t') d''^ , (2.6) 



where the Einstein relation (2.6) guarantees that the 
equilibrium probability density is given by 



/P[{V?}]e-^[W?>] 



(2.7) 



Following thej-dynamical perturbation theory devel- 
oped by JansserO and De Dominicisjlj we want to cal- 
culate the dynamical properties of our system, e.g. the 
dynamical correlation functions. 



Fir st, t he stochastic forces are eliminated, using equa- 



tion (2^) and the Gaussian distribution for the stochas- 
tic forces After a Gaussian transformation and the 
introduction of auxiliary Martin-Siggia-RosecJ fields V'" j 
the non-linearities occuring in the initial functional are 
reduced. A perturbation theory analogous to the static 
theory can now be implemented on the basis of the path- 
integral formulation. We define the generating functional 

Z[{/i"}, (X I V[{^i,^}]V[{^^}] 

X e"'[{C}:{C}]+/'i''a;Jdti:„(h° C-hh"C) ^ (2.8) 

where the resulting Janssen-De Dominicis functional J = 
^0 + Jint is split into the harmonic part Jq and the inter- 
action part Jint, 
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Ao ^o"(k,c^)C(-k,-c^) 



V'o"(k,c^) + Ao (ro + A:2) i^^(-k,-c^) 



(2.9) 



J^nt[{k}AK}] 



"An Uo 



X V ^o"(kl, ^l) ^o"(k2, C.2) (k3, C.3) V'f (k4, iUi) 



al3 



(2.10) 



The N-point Green functions G^^a^a(k, w) and cumu- 

lants C can be derived by appropriate derivatives of Z 
and InZ with respect to the sources and . Thus 
the standard scheme of perturbation theory can be ap- 
plied. Further details can be found in textbooks (Refs. 
and in Refs. 

In addition we want to list some important relations 
that will be useful for the discussion. The dynamical sus- 
ceptibility gives meaning to the auxiliary fields by noting 
that it can be represented as a correlation function be.- 
tween an auxiliary field and the order parameter fieldta 



X?''(x,t;x',t') 



5(C(x,t)) 



Sh0{x',t') 
= (Vo"(x,t)Ao ^f(x',0> 
Its Fourier transform 

X?'^(k,c.) = AoG,^„^,(k,c.) . 



(2.11) 

(2.12) 
The 



is associated with the Green functions G^^^^ff 
fluctuation-dissipation theorem relates the correlation 
function of the order parameter flelds and the imaginary 
part of the response functionlla 



Go^a^3(k, Cj) 



(2.13) 



which will enter the calculation of the NMR relaxation 
rate. E.g., considering only the harmonic part Jq of the 
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dynamical functional and carrying out the functional in- 
tegration gives 



2Ao 



[Ao(ro +P)]2+w2 



(2.14) 



Finally we want to introduce the vertex functions 
r„ 7a IB, which are related to the cumulants through a 

v*^' rir3 

Legendre transformation. For examplalJil3 



G" . ,„(k,t^) = - 



1 



O l/f*^ ip' 



,(-k,-a;) 



(2.15) 



The vertex functions are entering the explicit calculation 
of Z factors and the susceptibility in the renormaliza- 
tion group theory. The advantage of working with ver- 
tex functions is that they are represented by one-particle 
irreducible Feynman diagrams only. 



III. iVMii-EXPERIMENTS AND SPIN-LATTICE 
RELAXATION 



Quadrupolar perturbed nuclear magnetic resonance 
(NMR) can be used to study the dynamics i-of phase 
transitions from a TV to a IC modulated phaseJla In this 
method the interaction between the nuclear quadrupole 
moment Q and the electric-field gradient {EFG) V is the 
dominant perturbation Hq of the Zeeman Hamiltonian. 
Thus in the corresponding Hamiltonian 



Ti- = Hz + Ti-Q 



(3.1) 



next to the dominating Zeeman term Tiz one has to con- 
sider the quadrupole interaction Hq — ^ J2j k QjkVjk as 
a perturbation. The quadrupole moment operator Qjk 
is ccg^d linearly to the EFG tensor Vjk at the lattice 

The fluctuations of Vjk can be expressed via order param- 
eter fluctuations, because of the dor nin i an t linear coupling 
of the EFG to the order parameteii230 



5V,, (x, t) = Au-j [5^^{t) + i<5V"^(t)]' 



(3.2) 



with spin / = | in three dimensions are given byESS 



oc 



3x?'^(k,t^L) 



BZ <^L 



1 



Go^=^^(k,cjL)dfc , 
(3.3) 



with the spectral density of the EFG fluctuations 



V,,{m,*{t + T)e-''^^dT . (3.4) 



Measuring the spin-lattice relaxation thus yields infor- 
mation on the susceptibility of the local fluctuations of 
the order parameter. 

The spin-lattice relaxation was studied in great detail 
by -means of echo pulse methods both below and above 
r/.Ea'tj Below the critical temperature T/, it is possible 
for the prototypic system Rb2ZnCl4 to identify the re- 
laxation rates I/IY^ and l/T^, dominated in the plane- 
wave limit bji, the amplitudon and phason fluctuations, 
respectively.E3 Therefore, the dynamical properties of the 
order parameter fluctuations can be studied below the 
phase transition as well, and separately for the two dis- 
tinct excitations. 



IV. HIGH-TEMPERATURE PHASE 

In this section, the critical behavior of the incommen- 
surate phase transition above Tj will be investigated. On 
the basis of scaling arguments and the use of critical 
exponents, calculated within the renormalization-group 
theory for the XY model in three dimensions and model 
A, the temperature dependence of the NMR relaxation 
rate is analyzed in the first subsection. Next we study 
the crossover scenario of the temperature dependence of 
the relaxation rate in the second subsection by means of 
the renormalization group theory. Comparison with ex- 
perimental data is made, and we comment on the width 
of the critical region. 



A. Scaling laws for the relaxation time 



We now briefly sketch how the OP fluctuations deter- 
mine the relaxation rate. The spin-lattice relaxation de- 
scribes the return of the nuclear spin magnetization M in 
the direction of the external field back to its themal equi- 
librium value following a radio frequency pulse. l3 During 
that time the energy of the spin systenL-is transferred 
to single modes of the lattice fluctuations.Ej Because the 
EFG fluctuations can be written as OP fluctuations, the 
spin-lattice relaxation is determined by the spectral den- 
sity functions of the local OP fluctuations at the Larmor 
frequency lu = ujl- The transition probabilities for nuclei 



Above the phase transition, the thermodynamical av- 
erage of the order parameter components is zero. Because 
the structure of the correlation function of the order pa- 
ram etcr d oes not change dramatically above Tj (see sec- 



tion [VB ) and the calculation of the relaxation rate 1 /Ti 
involves the integration of the correlation function over 
all wavevectors, we will derive a form of the correlation 
function using scaling arguments. Thus we are able to 
discuss the universal features of the relaxation rate be- 
havior when approaching T/ from above. In the har- 
monic approximation we immediately get the correlation 
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function, which turns out to be the propagator of our 
functional Jq, [see Eq. ( ^.14| )] 



Go(fc,c^) 



2Ao 



[Ao(ro 



(4.1) 
(4.2) 



The suffix o will be omitted in the following discussion of 
this subsection, because no renormalization will be con- 
sidered here. 

We want to exploit our knowledge about the critical 
region. The static scaling hypothesis for the static re- 
sponse function states 



X{k) = ■ x{x) ■ k-^+'^ 



(4.3) 



with the scaling function x, a constant prefactor A^, the 
scaling variable x — (fc^)^^ denoting the correlation 
length), and the critical exponent r]. Neglecting a fre- 
quency dependence for the kinetic coefficient (Lorentzian 
approximation), the dynamic scaling hypothesis for the 
characteristic frequency of the OP dynamics states 



u^{k) = A(fc)/x(fc) ^ fc" , 
and we can deduce 

\{k) = Ax ■ \{x) ■ ^-2+" , 



(4.4) 



(4.5) 



with X{x) being the scaling function for the kinetic coef- 
ficient and A\ a c onst ant prefaetor. Notice that for fixed 
wavevector k Eq. 4.4 leads toE3 



(4.6) 



The correlation function G(fc, lu) can now be rewritten in 
scaling form 



G{k,u) 



A. 



1 



kz+2-n 



(4.7) 



with 



A = — r^- 

Ax ■ 

/(d), a;) = x{x) ■ 



Ax k^' 
X{x)/x{x) 



1 

k~e 



[X{x)/x{xW+uj' 



(4.8) 



where the Lorentzian line shape is retained. Above T/, 
this not a very crucial approximation, because the shape 
of the correlation function does not change in a first-order 
renormalization group analysis, as we will see in the next 
section. To calculate the relax ation rate 1 /Ti , one has to 
evaluate the integral [see Eq. (|3.3|)] 



1 

7^ Of 
Ti JO 

cx A • 



1 



k'^ G{k, ijj[^)dk 



(4.9) 



k^dk k-'--'^+'^f{ku:l^'\kS) . 



BZ 



— 1/z 

With u = kujj^ and v = k^ we introduce new variables 



y/vF^rV^ ~ k\J ^j^^^ + 



(4.10) 



and 



taniy9 



(4.11) 



This leads to the relation 

' " I,[fig,^)) , (4.12) 

where the integral Ig does not contribute to the leading 
temperature dependence. 

The temperature dependence of the relaxation rate can 
now easily be found in the limits where the Larmor fre- 
quency or the frequency, of the critical fluctuations, re- 
spectively, dominate the integral and its prefactor. 



1. Fast-motion limit (uiL/i^f^ 



For temperatures very far above the critical temper- 
ature T/, the characteristic frequency is larger than the 
Larmor frequency. Thus the temperature dependence of 
the OP fluctuations determines the temperature depen- 
dence of the relaxation rate; the ialue of the Larmor 
frequency should not play any role.tJ For the integral 



1 



A 



^ iz-i-n)/2 



(4.13) 



we obtain with tanip = const, [see Eq. 4.11 



T -Ti 
T 



~IJ-{Z-1-T]) 



(4.14) 



Taking the values for the critical exponents from Table 
0, we find 



Ti 



T-Ti 



-0.663 



(4.15) 



This can be compared with_the experimental results 
for Rb2ZnCl4 by Holzer et al.EHl, who found for the lead- 
ing scaling behavior of the relaxation rate, following the 
temperature-independent region, the exponent —0.625. 
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2. Slow-motion limit (ljl/i-Oip ^ i) 

In the vicinity of the critical temperature Tj critical 
slowing down will occur. This means that the character- 
istic frequency is approaching zeco and will fall below 
the value of the Larmor frequency.^ Thus, the charac- 
teristic time scale of the OP fluctuations is slower than 
the experimental time scale. For the temperature depen- 
dence of the relaxation rate 



1 



rocA.(e^ [(r7-L)'/^ + i]) 



^ (z-l->))/2 



In 



we now obtain 



Tl 



■ const 



(4.16) 



Taking the values for the critical exponents from Table || 
again 



-0.49 



(4.17) 



This., is in good agreement with the experimental 
resulio for Rb2ZnCl4 that the value of the relaxation 
rate near T/ scales as ujj^'^'^ for different Larmor frequen- 
cies. 

We want to stress that the transition from the fast- to 
the slow-motion limit is a property of the integral enter- 
ing the calculation of the relaxation rate. Because the 
susceptibility is evaluated at fixed ujl, there exists for 
a lower boundary for the integral T Tj. This means 
that the transition from the temperature-dependent to 
the temperature independent behavior near Tj is fixed 
by the scale lo^. It should also be mentioned that our 
results reproduce those obtained earlier by Holzer et al. 
(see Ref. pO| ), when the van Hove approximation for z is 
used (z « 2). 



B. Renormalization-group analysis above Ti 

Our investigations concerning the critical behavior 
above Tj in the last section led to fair agreement with 
experimental data. It was possible to gain the critical ex- 
ponents for the frequency and temperature dependence 
of the relaxation rate in a quantitatively accurate way. 
Furthermore, we obtained a qualitative understanding of 
the transition of the relaxation rate 1/Ti from the slow- 
to the fast-motion limit. This transition is caused by the 
characteristic frequency of the order-parameter dynamics 
uj^ approaching zero, i.e., the critical slowing down near 
T/. This renders the Larmor frequency lo^ the dominat- 
ing time scale, 1/Ti becomes temperature-independent 
near Tj. 

At this point we want to consider what happens upon 
leaving the region near the transition temperature and 
going to higher temperatures. Our results (4.13) for the 



fast-motion limit are based on the assumption that fluc- 
tuations are very important. We used non-classical expo- 
nents and scaling arguments valid in the critical region. 
The question arises how large this region of temperature, 
where the system displays the non-classical behavior cal- 
culated in the last section, will be. Increasing tempera- 
ture diminishes the effect of fluctuations. One would ex- 
pect that at some temperature Gaussian behavior should 
emerge. We shall apply the renormalization group theory 
to one-loop approximation in order to describe the tran- 
sition from the fluctuation-dominated behavior near Tj 
to a temperature region where the mean-field description 
should be valid. It is obvious that the properties of the 
integral, responsible for the transition between the slow 
and the fast motion limit, will account for the crossover 
of the leading scaling behavior when the temperature is 
increased. In order to discuss the crossover within this 
analysis a modified minimal-subtraction prescription is 
employed. ThiS|-scheme was first introduced by Amit 
and GoldschmidtEj and subsequently explored by Lawrie 
for the study of Goldstone singularities .Lj It can com- 
prise exact statements in a certain limit. Below Tj, this 
is the regime dominated by the Goldstone modes alone; 
in the region above T/, which we will consider in this 
section, the mean-field result is used. In this scheme, in 
addition the standard field-theoretical formulation of the 
renormalization group is neatly reprod-uccd. Following 
the arguments of Schloms and DohmE3 and Ref. ^ we 
can refrain from the e expansion, with 



(4.18) 



defining the difference from the upper critical dimension 
of the model. This is motivated by the following. 
Above Tj the Gaussian or zero-loop theory becomes exact 
in the high- temperature limit. The critical fixed point, 
i.e., the Heisenberg fixed point, dominating the behavior 
of the system near the critical temperature is calculated 
to one-loop order. The main interest here lies in the 
crossover behavior between these two fixed points, which 
is calculated to one-loop order, too. Thus no further ap- 
proximations are necessary to be consistent. 

Very close to the critical teaaperature Tj an e- 
expansion or Borel resummationE3 would of course be 
inevitable in order to obtain better values for the critical 
exponents. A description of the generalized minimal sub- 
traction scheme is for example given in Refs. |35| , [l5| , |l6| and 
|37| . A crossover into an asymptotically Gaussian theory 
is described by this method in Ref. |l^. 



1. Flow equations 

Our aim is to calculate the wavenumber and frequency 
dependence of the susceptibility to one-loop order. The 
field renormalization is zero to one loop order. Thus we 
will not take into account corrections to the static ex- 
ponent rj and corrections to the mean-field value of the 
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dynamic exponent z « 2. This leajies To and Uo as the 
only quantities to be renormalized.cJ 

There is a shift of the critical temperature from the 
mean-field result To/ to the "true" transition tempera- 
ture T/. In order to take this shift into account, a trans- 
formation to a new temperature variable, being zero at 
the critical temperature T/, is performed. This new vari- 
able will be denoted again as Tq. The only renormalized 
quantities are then written as 

-2 



Here, the geometric factor Ad is chosen0 

r(3 - d/2) 



as 



Ah 



(4.19) 
(4.20) 

(4.21) 



For the non-trivial Z factors one finds in the general 
ized minimal subtraction procedure (see App. ^) 

" + 8 _ 1 

-UoAdfi 



Zu^l 



1 



6e 



(1 



-UoAdH" 



1 



(4.22) 
(4.23) 



6e """"^^ (l + ro/M2)^/2 

Setting To = the familiar renormalization constants 
for the n-component 0'* model are recovered. In gen- 
eral here, however, the Z factors are functions of both 
Mo and To £3 In the next step the fact that the un- 
renormalized A^-point functions do not depend on the 
scale /i is exploited and the Callan-Symanzik equations 
are derived.EEl The idea behind that is to connect via 
the renormalization-group equations the uncritical the- 
ory, which can be treated perturbationally, with the crit- 
ical theory displaying infrared divergences. The result- 
ing partial differential equations can be solved with the 
method of characteristics (/i(0 = mO- With the defini- 
tion of Wilson's flow functions 

Cr(0 = /i:^ 



dpi 
d_ 
dp 



In — 





(4.24) 



(4.25) 



we procee d to the flow d ependent couplings t[1) and u{l) 
(see Eqs. and |^) 



Bu{l) 



dl 



(4.26) 
(4.27) 



given by the flrst order ordinary differential equations 



drjl) 
dl 



Til) 



-U[l) 



1 



-U{1) 



[1 + t(0]1+^/2 



[1 + t(0]1+^/2 



(4.28) 



(4.29) 



and the initial conditions t(1) = t and u(l) = u. The 
asymptotic behavior is determined by zeros of the (3 func- 
tion, giving the flxed points of the renormalization group. 
Here, we find the Gaussian fixed point u*q — with 

= — 2 and the Heisenberg fixed point — 
with C,^ "rT^ + £■ These fixed points are of course 
well-knownpj but in the generalized minimal subtrac- 
tion scheme it is now possible to describe the crossover 
between these two fixed points. We are interested in the 
theory in three dimensions and will henceforth discuss 
this case (e = 1). 

First we investigate the crossover of the t(/) fiow. It is 
possible to recover the universal crossover in the flow by 
plotting t{1) against the scaling variable (compare Refs. 



I 



T(l) 



1/(2- 



(4.30) 



In Fig. |i| the effective exponent for the a;-dependence 
of t{1) is depicted for ten different values of t(1) [with 
flxed u{l) and n = 2], coinciding perfectly. There is a 
crossover from the region / with the exponent —2 to 
the region I 1 with the exponent — 2 + (n -I- 2)/(n -I- 8). 

Next we find, with the scaling variable x oc {k£,)~^, the 
effective exponent VcS of the temperature dependence of 
the correlation length 



t(/) cx r^^""'' 




n+2 
ri+8 







(4.31) 



Thus, with the generalized minimal-subtraction scheme 
we can describe the crossover from the non-classical criti- 
cal behavior to the Gaussian behavior, e.g., as a function 
of the temperature variable t. 



2. Matching 



To one-loop order, only the tadpole graph enters in 
the two-point function (see r^^^^^ in App. |^) shifting the 
critical temperature as stated above. Thus the suscep- 
tibility does not chang e its form a nd t he renormalized 
version reads with Eq. ( |2.12| ), Eq. ( |2.15D and App. ^ to 
one-loop order. 



(4.32) 



Yet what we gained in the last subsection was the tem- 
perature dependence of the coupling constants. We have 
to take into consideration this temperature dependence 
in order to discuss the changes resulting from the fiuctu- 
ation corrections. One now has to ask the question: how 
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does the flow enter the physical quantities measured in 
an experiment? 

With the flow dependence of the couphng constants, 
t he relaxation rate to one-loop order becomes [see Eq. 

(P)] 



1 OC / fc2^fcSXi.(k,C.,) 

Ti J ujl 
= ( k^dk 



1 



^l/X^ + [k^ + T{l)^l^^Pf A 



1 1 

III A 



k^dk- 



[fc2 



(4.33) 



where k = k/ and uj^ = ujl/ X^'^P . 

Keeping r{l) fixed (it is set to 1) the relaxation rate 
1/Ti is proportional to for large I. When I approaches 
zero a constant value of the integral and hence of the re- 
laxation rate 1/Ti will be reached, because of the fixed 
time scale 1 /(D^ . The physical reason is that in the slow 
motion limit the characteristic time scale (l/wi^) becomes 
larger than the experimental time scale In Fig. |^ 

the logarithmic I dependence of the integral 



h = 



aiogTi 

dlogl 



9 log I — — / k'^dk — 



i + (P + l)2 



dlogl 



is plotted against the scaling variable 
X = I ■ 



(4.34) 



(4.35) 



We regain the transition from the ^-independent regime 
(/ — » 0) (and therefore of the temperature-independent 
regime) to the regime Ii (x I {I ^ oo). This corresponds 
to the transition from the slow-motion to the fast-motion 
limit; a change of I is aquivalent to a change of Lu^p. 

However, we are rather interested in the dependence 
of the relaxation rate from the physical temperature r(l) 
than from the flow parameter /. This may be obtained as 
follows. Knowing the solution of the flow equations t{1) 
and u{l), we can flnd a h for a given t(1) that fulfills 
the equation t(/i) = 1. Inverting this relation, r(l) for 
a given li with t{Ii) = 1 can be found. It is not possible 
to write down an analytical expression, but numerically 
this relationship is readily obtained. Thus we are led to 
l/ri(r(l)) = l/ri(Zi[r(l)]). To connect the theory in 
a region where the perturbation expansion is valid with 
the interesting region, we match the temperature vari- 
able r{l) to 1, thus imposing the crossover behavior of 
the flow r{l) to the effective exponent of the relaxation 
time 1/Ti. 

In Fig. ||the resulting Ti(T) dependence is used to fit 
experimental NMR data for Rb2ZnCl4 (measured points 
are indicated by circles) from Mischo et al.ES Two pa- 
rameters have to be fixed in the theory. First, the pref- 
actor relating the relaxation rate and the integral over 



the imaginary part of the susceptibility in Eq. (4.33) 
must be determined. Thus, the value of l/Ti{T — Tj) is 
set. The second parameter is the scale of ujl compared 
to the coupling A/i. With this the relative temperature 
AT, where the transition between slow and fast motion 
limit takes place, is adjusted. 

The two fit curves presented in Fig. |3|show a crossover 
to the mean-field regime starting at AT « lOi^T (a) and 
at AT w 5K (b). A very good agreement, not only for 
the transition from the slow to the fast motion limit, but 
also for the high-temperature behavior is found in the 
second case. We want to discuss this issue now in more 
detail. 



3. Width of the Critical Region 

Some experiments report large regions in which non- 
classical exponents for the temperature dependence of 
the relaxation rate are observed. E.g., in Ref. ^ the 
range above Tj where non-classical exponents are found 
is AT « lOOK. These findings haye-ta be understood by 
means of the Ginzburg-LevanyukOcJ argument, which 
states that only near to the critical temperature the non- 
classical critical exponents should be valid. Fluctuations 
should contribute only near the critical point and change 
the mean-field picture there. 

The property of the integral quantity 1/Ti, in the re- 
gion where a crossover between the non-classical criti- 
cal exponents and the mean-field exponents occurs, was 
studied in the last subsection. We now comment on the 
four regions that can be identified and are indicated by 
numbers 1 ... 4 in Fig. |[ 

Very close to T/, there is a temperature-independent 
region (1), because of the dominating scale lul- Here, 
the probing frequency lol is too fast to grasp the criti- 
cal behavior. Upon going to higher temperatures, after 
a transition region (2), a temperature dependence with 
non-classical critical exponents emerges (3). For even 
higher temperatures one finds a crossover to the mean- 
field exponents, in regime (4). In Fig. I this crossover 
takes place between AT w 5K and AT « 20K. From 
Fig. H, we find that the crossover at lower tempera- 
tures, here starting at AT « 5K, leads to a better fit of 
the experimental data. Thus the reported large reeipn, 
where supposedly non-classical exponents are found,E3 is 
in our opinion not an indispensable conclusion that can 
be drawn from the experimental data. The plausible sce- 
nario of an extended crossover regime beyond the truly 
asymptotic region of width AT w 5K is in fact in perfect 
agreement with the data. 

As this is not a universal feature other scenarios are 
possible. It may happen that the scale of is very 
large and thus only the Gaussian exponents are found. 
We omitted the contribution p£ higher Raman processes, 
as discussed by Holzer et al.,E2l leading to an additional 
T^ dependence for the relaxation rate. These would bend 
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the curves downward even more and explain the devia- 
tion present at the highest measured temperatures. Not 
taking these additional contributions into consideration, 
however, clarifies the crossover aspect. 



V. LOW-TEMPERATURE PHASE 



This section is devoted to the incommensurate ordered 
phase below the critical temperature Tj. In the 0{n)- 
symmetrical model a spontaneous breaking of a global 
continuous symmetry occurs and the expectation value of 
the order parameter becomes nonzero. Now parallel and 
perpendicular fluctuations with respect to the nonzero 
order parameter have to be distinguished. As a con- 
sequence there appear n — 1 massless Goldstone modes 
which lead to infrared singularities for alLtemperatures 
below Ti in certain correlation functionsJIlu We investi- 
gate how these Goldstone modes influence the dynami- 
cal properties of the quantities we are interested in, e.g., 
the NMR relaxation rate. To do so, we first derive the 
dynamical functional appropriate below Tj. In the fol- 
lowing subsection some comments about the Goldstone 
anomalies are made. We will then treat the dynamics of 
the fluctuations parallel (amplitudons) and perpendicu- 
lar (phasons) to the order parameter. Again a renormal- 
ization group calculation to one-loop order is presented. 
We will discuss the dynamical susceptibility before eval- 
uating the integrals leading to the relaxation rate. In the 
last section, we compare with experimental data. We 
also comment on the existence of a phason gap. 



A. Dynamical functional 

Let us assume that the spontaneous symmetry break- 
ing below T/ appears in the 71th direction of the order pa- 
rameter space. As usuai, new fields 7ro,Q; = l,...,n— 1, 
and (To are introducedlij 



71-0 _ 
(To + 4>o 



(5.1) 



(5.2) 



(5.3) 



Thus (To corresponds to the longitudinal, and tt" to the 
transverse fluctuations. 

Ins erting these transformations into the functional 
(2.10) leads to a new functional of the form J = Jq + 



with 

«) - ((To) - . 

The order parameter is parameterized as 



J'int + Ji + const withQ 

^Ao <(k,C^)7f„"(-k,-Cc>) 

'- Q. 

+ Ao (To(k, cj) ao(-k, -uj) 
2 

^<(k,c^) [loj + K (^° + ^ + ^')] <(-k,-u;) 

a 

r / 3 ihj^ \ 1 
(To(k, w) ibJ + \o [toA 2^^^^) o'o(— k, — Cj) , 

(5.4) 



Jint[{TTo},do,{T^o},CFo 



\ Ao U. 

6 



^ TTo (kl , Wl ) TTo (ka , CJ2 ) Trf (kg , W3 ) 7rf (k4 , CJ4 ) 

■ a/3 

- ^ ^o (kl , CJl ) TTo (k2 , W2 ) (To (ks , ) (To (k4 , CJ4 ) 

a 

- ^ (To (ki , ) TTo (k2 , W2 ) TTo (ks , ) CTo (k4 , W4 ) 

a 

- CTo (ki , ) (To (k2 , ^2 ) (To (ka , ) CTo (k4 , W4 ) 

J2 2 <(kl, C^l) <(k2, C^2) (To(k3, ^3) 

a 

- ^ (Jo (ki , Wi ) Vr" (k2 , W2 ) Tt" (kg , W3 ) 

a 

- 3 ao (ki , wi ) CTo (k2 , 0^2) (To (k3 , uj^) 



(5.5) 



and 



T r~ 1 X /3 f ^o 

(To = - Ao 1/ Too \ro + — 

\ Uo \ 2 



(To(— k, — a))(5(k, L 



(5.6) 



Equation ( |5.2| ) (((To) — 0) yields a perturbative identity 
that gives thCi-Kelation between Vo and TOo, reading to 
one-loop ordeiO 



n — 1 
6 



1 



kro + ^ + k^ 



2~"./fe ro + ^ + k^ 



(5.7) 



In the following To is replaced by TOo. Notice that by 
using the varia ble TO q , the shift of Tj is already incorpo- 
rated [see Eq. (|5.3|)]. We can now write down the basic 
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ingredients needed to apply the recipe for the dynami- 
cal perturbation theory below T/. The emerging prop- 
agators, vertices, and counterterms are listed with their 
graphical representation in Figs. |[ ^ and (see Ref. 16). 



B. Goldstone theorem and coexistence limit 

As mentioned before, the particularity of the 0{n)- 
symmetric functionals below the critical temperature is 
the occurrence of |Goldstone modes in the entire low tem- 
perature phascElQ Because no free energy is required for 
an infinitesimal quasistatic rotation of the order parame- 
ter, the transverse correlation length diverges in the limit 
of zero external field. The corresponding massless modes 
are the Goldstone modes,cl in this context called pha- 
sons. They are manifest in non-analytical behavior of 
correlation functions, for example the longitudinal static 
susceptibility diverges and chang^^|ts..lfiading behavior 



from being proportional to fc t 



Xl 



-i(k,0) cx k" 



(5.8) 



Before discussing the details of the renormalization the- 
ory below Tf, we summarize some important aspects, 
which explain why below Tj an e-expansion can be 
avoided. For more details see Ref. ^ 

Leaving the critical temperature region T w T/, which 
is characterized in the non-perturbed case by TOq = 0, 
and lowering the temperature, means that the fluctua- 
tions of the longitudinal modes (amplitudons) become 
negligible, because these modes are massive (mo). In 
contrast the phasons remain massless and hence their 
fluctuations will dominate. 

Yet a different way of describing the dominance of the 
fluctuations of the Goldstone modes is to consider the 
spherical model limita n —> oo. In this case of "maximal" 
symmetry breaking, the Goldstone modes are weighted 
with the factor rt — 1 — )■ oo. 

As mentioned above, below Tj coexistence anomalies 
are present. They arise from the fact that the n—1 trans- 
verse modes are massless. In the limit fc ^ and w — > 
for T < Tj this manifests itself in typical infrared diver- 
gences. An important point to remember is that in order 
to gain these coexistence anomalies, one can also study 
the case mo — > cx). In the renormalization scheme it is 
shown that the flow of the mass parameter nio tends to 
infinity as the momentum and frequency tend [to zero. 

From this it is plausible, and was also provedcJ, that in 
the coexistence limit the result for the two-point vertex 
functions are identical with the results arising from the 
spherical model limit n — > oo. 

These findings render an e expansion unnecessary in 
the coexistence limit (fc ^ 0, w ^ at mo > 0), because 
the asymptotic theory (the spherical model) is exactly 
treatable and reduces to the zero- and one-loop contri- 
butions. Of course, one has to make sure that the prop- 
erties of the asymptotic functional will be reproduced in 



the respective limit. Within the generalized minimal sub- 
traction scheme this isj)ossible. As stated in subsection 
IV B l La wrie's methodts and its dynamical extension in 
Ref. |16| lead beyond these limits and allow for a detailed 
study of the crossover behavior. 

The behavior of the correlation functions in our case 
is driven by the crossover between the three fixed points 
present below Tj. Besides the Gaussian fixed point one 
finds the Heisenberg fixed point 



[UH = 6e/(n + 8)] 
and the coexistence fixed pointEl 



(5.9) 



(5.10) 



We will again employ the generalized minimal subtrac- 
tion scheme to study the crossover between these fixed 
points. 



C. Renormalization group analysis below Tj 

1. Flow equations 

Below Tj, using only the one-loop diagrams, again 
the field renormalization vanishes. Hence, the only non- 
trivial Z factors are the ones for the temperature scale 
and the coupling constant: 



2 2 -2 

m = Z„ m„/i 



(5.11) 
(5.12) 



Because we use m instead of r an important relationship 
can be-stated, which is true independently of the loop 
order,Ei 

Zrn • Zcr = Z^ ■ (5.13) 

To one-loop order {Z„ — 1) we find (see App. ^)lll 

Zu = Zjn =1 + ^ UoAdfi~'' 

6e 

Here, the contribution of the transverse loops lead to dif- 
ferent divergences manifest in the cha nge of the Z factors 
compared to those above T/ [see Eq. ( 4.23 )] . We recover 
the familiar renormalization constant in the critical re- 
gion by setting mo = 0. When considering the coexis- 
tence limit mo oo, the weight of the effective critical 
fiuctuations is reduced from n + 8 to n — 1, the number 
of Goldstone modes. 

Asymptotically (m —^ oo) the Z factors are exact. In 
the crossover region they—are an approximation to the 
order u^/ {1 + rrio/ ^l^y/^lS 
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From this we directly derive the flow-dependent cou- 
plings 



dl 



dujl) 



u{l) 



2 [l + m(02]i+^/2 



+ 



u{l) 



2[l-Fm(/)2]i+^/2 



Wilson's flow equations /3„ and Cm now read 



dl 



(5.15) 

(5.16) 

(5.17) 
(5.18) 



Three flxed . pnipts have now to be taken into 
consideration.E£lE3 Next to the Gaussian flxed point 
u*Q = with = —2, we find in the critical limit 

(mo 0) the infrared-stable Heisenberg fixed point 
u*u = &£/{n + 8) with C,^fj = -2 -I- e. In the coexis- 
tence limit mo — *■ cx3, we find in addition to the still 
ultraviolet-stable Gaussian fixed-|point the coexistence 
fixed point, identified by Lawriejl3 = 6e/(n— 1) with 
QnH = ^2 + e, which is infrared-stable. Thus m{l)^ 
diverges asymptotically for Z — > as if e < 2. In- 

deed, the coexistence limit is described by a divergent 
mass parameter. 

In Figs. H and |^ the fiow for m{l) and u{l) is plotted. 
We find for the flow u{l) a crossover between the coexis- 
tence flxed point, inversely proportional to the number of 
Goldstone modes (n — 1), and the Heisenberg flxed point 



uil) 



6e/{n~l) I 
6e/(n + 8) / 



> 
1 . 



(5.19) 



That means that for m(l) <C 1 the coexistence limit is 
not approached directly for I 0, but for a while the 
flow stays near the Heisenber g iyie d point regime. The 
scaling variable for u{l) is heret£Hl3 



I 



m(l)2/(2-£) 



(5.20) 



again leading to perfectly coinciding curves when plotted 
vs. X. 



From the relation stated in Eq. ( ^,3 ) one can deduce 

riantlj 



the renormalization-group invariantt 



Mlfp-e ^ "l(l)2 



uil) 



U{1) 



(5.21) 



which immediately gives us the scaling of m(Z) that can 
be observed in Fig. ra. 



to(/)2 cx r^/"-" with = 



1 



2-e 



n+2 
n+8 







(5.22) 



Notice that the value of l/i'cff in the first case is the same 
as l/v for the spherical model. The mass parameter m 
diverges for I 0, m(l)'^ cx /"^^^ with e < 2. 

From now on we will concentrate on three dimensions 
(e — 1) and n = 2. 



2. Matching 

In order to discuss the static susceptibility we use the 
matching condition jil = k. This relation connects the 
dependence of the renormalized quantities on the mo- 
mentum scale fi with the k dependence, in which we are 
interested. 

For the integrated value, we are interested in the tem- 
perature dependence rather than the dependence on the 
flow parameter I or the wavevector k. Thus, after inte- 
gration we have to again match the resulting Z-dependent 
relaxation rate with the physical temperature (see section 
fVBl). 



D. Susceptibility 

In order to determine the renormalized dynamical sus- 
ceptibility, we evaluate the one-loop diagrams for the 
two-point cumulants, which can be easily derived—from 
the one-loop vertex functions listed in Appendix Alll3 Be- 
low Tf, the structure of the susceptibility does change to 
one loop order, compared to the mean-field results. We 
write in a general form 



x;i/i|(k,c.) = — :i + fc2 + /o^/"(k,c.) 



(5.23) 



with the self-energy /o*"^" containing the contributions 
of the one- loop diagrams. The explicit form of /c^^" is 
gained from the calculation of the two-point vertex func- 
tion in App. ^ and Eqs. ( p.l2D and ( pTS] ). We then ob- 
tain the renormalized susceptibility X\\' ± by inserting the 
Z factors with the fiow dependent coupling constants u{l) 
and m(l). Because no field renormalization is present, we 
can replace Ao with A to this order. This is because via 
the fluctuation-dissipation theorem ( 2.13| ), Zx and the 
held renormalizations are connected. 
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1. Amplitudon modes 

In d — 3 the longitudinal susceptibility characterizing the amplitudon modes is given bj^l 

1 u(l)m{lf 



n — \ { TT . tuj/X \ 

— - — - + arcsm — — — +3 arcsm ■ 



iLu/Xfjr I 



2 12 



{[k^ ^iuj/X]/fi^ Pf +4m{l)^k^/fi^f 



1/2 



arcsm ■ 



(fc2 -Ztj/A)/^2;2 



([fc2 - iu/X]/fi^ Pf+A m{l)^ fc2 /^2 



1/2 



(5.24) 



First we want to discuss some limits in order to become 
acquainted with this complex form of the susceptibility. 

It is important to notice the change of the structure 
of the RG susceptibility that results from the one-loop 
contribution of the perturbation theory. To clarify this, 
we state the asymptotic susceptibility, which is evaluated 
for non-zero frequency (k ^ 0, w > 0) 



X|| 



(5.25) 



fc2 _itj/A + m(l)2 . k 



1 



■9{k) 



with a constant a and a function g{k) that is regular for 
g — > 0. From this limit, it becomes clear that we have to 
expect changes of the scaling behavior. 

In the coexistence limit (a; = 0, k 0) we recover the 
exact asymptotic result (d = 3, e = 1) 



Xf oc 



(5.26) 



displaying the coexistence anomaly. When keeping the 
frequency lu > fixed (k — s- 0, w 0) the imaginary part 
of the susceptibility approaches a constant value 



h{Lu) 



(5.27) 



where h{u!) is a function of lu only. We can now turn 
to the full susceptibility. The imaginary part of the sus- 
ceptibility is plotted for different temperatures in Fig. 
p^ . The structure of '^xf changes dramatically as com- 
pared to the mean-field result, as to be expected. The 
contributions of the phason and amplitudon loops are 
given by the terms in brackets of Eq. (5.24). They give 
rise to a qualitatively different behavior of xf - Different 
scaling regions can be identified. Expanding the imagi- 
nary part of xf yields analytical expressions for the scal- 
ing regions, as listed in Table ||. While the k 00 
and k —^ behavior reproduces the mean- field result, 
the correct treatment of the Goldstone anomalies lead 
to an additional k~^ behavior in the intermediate region 
y/uj/X < k < m(l). A plateau appears for smaller k and 
temperatures far away from the critical temperature T/. 
The effective exponent k of the k dependence of "^xf is 
plotted in Fig. |ll[ One can therefrom easily identify the 
scaling regions presented in Table ||. 

The influence of the Goldstone modes is therefore to 
alter the k dependence of the susceptibility not only in 
the coexistence limit, but also in intermediate regions. 
In order to derive the temperature dependence of xji*', in 
addition the flows of m{l) and u{l) need to be considered. 



2. Phason modes 

For the transverse susceptibility characterizing the phason modes one findslii 



xl{K^) = 



1 



1 



u{l)m{l)/6 



k^-iuj/X\ (fc2-ia;/A)//x2p 



m(l) I TT 
2 — — - — - — — arcsin 
k/iil \ 2 



arcsin - 



-I- arcsin ■ 



iuj/XiJi^ P + m{iy 



((fc2 - iulX)l[L'^ P - m{Vff + 4 m{Vf k^ / p] ^'"^ 

~-iuj)/XiPP ~m{lf ^ 

((fc2 - iujlX)ly? P - m{lff + 4 m{lf k^ / ^jfl 



(5.28) 
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Here the problem lies in the cancellation of terms with 
respect to their k dep endence, hidden in the complex 
structure of Eq. ( 5.2§| ). Hence we start again with con- 
sidering the coexistence limit (k ^ 0, ^ 0): 



(5.29) 



which is easily found. The results for k ^ 0, w ^ are 
more difficult to obtain, because the arcsin-terms cancel 
their ^-dependence against each other. In two limits this 
can be done analytically. For m ^ one gets 



1 



Xj 



XoJ 



fc2 — iuj/ X 



(5.30) 



reproducing the mean-field susceptibility for the massless 
transverse modes. 

For m —> oo the arcsin-terms read as ^ ' ^ 
leading to 



-C2 



1 



Xj 



(5.31) 



6 



Here ci and C2 are constants. Thus in the two extreme 
limits, the temperature dependence vanishes, and only 
in between, for m{l) cx 0(1) can we expect a slightly 
temperature-dependent behavior. 

In Figs. ^ and |l^ the imaginary part of the trans- 
verse susceptibility and its effective exponent with re- 
spect to k are plotted for different temperatures. Notice 
that leaving the critical temperature leads to a temper- 
ature dependence. This is caused by the coupling of the 
amplitudon and phason modes. Yet as the temperature is 
further reduced, the temperature dependence disappears 
again. 



E. Relaxation rate 

In this subsection we study consequences for the re- 
laxation arising from the Goldstone anomal ies p resent in 
the susceptibility. As mentioned in Sec. Ill, in order 



to gain the relaxation rate we have to integrate over the 
imaginary part of the susceptibility. 

Because the transverse susceptibility is temperature 
dependent, also the relaxation rate, connected with the 
phasons, will be temperature dependent. This is of 
course not the case in the mean-field analysis. As dis- 
cussed in the last section, for T ^ Tj the susceptibility 
approaches the mean-field result, and thus the relaxation 
rate at the critical temperature is unaltered. 

For the relaxation rate connected with the amplitudons 
the changes are more subtle. Therefore, we collect all 
the contributions from the one-loop diagrams in a func- 
tion /(k, w), which can be interpreted as a k- and lu- 
dependent dimensionless self-energy. The susceptibility 
has now the following structure 



2;2 



(xf)"' =fc'-WA + /(k,u;)M"? 



(5.32) 



The dependence of / on fc and to is plotted in Fig. |TJ. We 
see that the real part of the effective mass / is decreas- 
ing for fc ^ 0. Thus, the Goldstone anomalies lead to a 
reduction of the real part of /. In the coexistence limit, 
uj ^ and s mall fc , the real part of / tends to linearly 
and relation (5.26) is recovered. The imaginary part is 
only fc-dependent for very small k. That means, when we 
integrate over all k the infiuence of the Goldstone anoma- 
lies can be interpreted as follows. The effective Larmor 
frequency is raised and the mass is lowered for small k as 
compared to the mean-field description. We can easily 
derive this easily from the longitudinal relaxation rate 
with f(k,Lu) taken into consideration 



3xf(k,c^L) 



k^dk 



k^dk 



1 



u) - 3/ll(k,<I>) 



p2 _ CJ/II (k, Cj)f + [fc2 + 3fi/ll (k, Cj)f 



1 1 

/iZ A 



k^dk 



1-S5/Il(k,ci)/Ci 



3/ll(k,c:;)]2 + [A:2 + 3f?/ll(k,c:;)]2 



(5.33) 



where again k — k/^l and Co — lol/X^'^P. When we 
compare this result to the mean-field result, the inter- 
pretation given above becomes clear. The relaxation rate 
is raised through the infiuence of the Goldstone modes. 
Both the transverse and longitudinal relaxation times are 
plotted in Fig. ^ Again we have compared our find- 
ings with experimental data, taken from Ref. In the 
low-temperature phase we have less freedom of choice in 
our theory, as the scale Ti(T = Ti) and the parameters 
are already fixed by their high-temperature values. Thus 
only one parameter is left to be adjusted. 

In the vicinity of Tj we find a temperature-independent 
region, because both the transverse and the longitudinal 
susceptibility become temperature-independent and ap- 
proach their mean-field values. 

The transverse relaxation time shows a slight tempera- 
ture dependence for temperatures further away from T/. 
If we use the identical choice of parameters as for the 
high-temperature phase, we find good agreement in the 
low-tcmperaturc phase as well. The temperature where 
the maximum value of the transverse relaxation time in 
our theory is reached is identified with the correspond- 
ing temperature in the experiment. This temperature 
dependence is due to the coupling between the phason 
and amplitudon modes. We want to emphasize that, 
in agreement-ipdth the analysis of ultrasonic attenuation 
experimentsjjll3 no phason gap has to be introduced to 
explain the experimental data for Rb2ZnCl4. However, 
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it is important to treat the influence of the Goldstone 
modes beyond the mean-field approach. 

For the longitudinal relaxation time, the crossover 
temperature represents additional an important scale. 
We used the same range AT « 5K as in the high- 
temperature phase for the plot in Fig. ^ Again good 
agreement between experiment and theory is observed. 
Both theory and experiment show two scaling regions, 
one above AT « 5K and one below. The qualitative 
behavior is correctly reproduced, but the quantitative 
agreement for the longitudinal relaxation rate is not as 
good as compared to the high-temperature phase. A pos- 
sible reason may be the following. We calculated the 
coupling of the transverse and longitudinal order param- 
eter fields to one- loop order. Below Tj, the coupling of 
the order parameter is changing the susceptibility in its 
structure, whereas above Tj nothing dramatic happens. 
One has to expect that below T/ this is of course only 
the first step beyond mean-field theory and the two-loop 
corrections might lead to quantitative modifications in 
the crossover region. Comparing the calculated trans- 
verse and the longitudinal relaxation rates below Tj with 
the experimental data is in agreement with this. A slight 
temperature dependence is not as sensitive as the scaling 
behavior of the longitudinal relaxation time, which has 
again two regimes due to the crossover scenario. 

As the characteristic features are reproduced correctly, 
we may say that we understand the complex temperature 
dependence below T/ in the context of the coupled order 
parameter modes and a careful treatment of the Gold- 
stone modes. Upon introducing the k- and w-dependent 
self energy / we could in addition provide a physical in- 
terpretation of the changes of the longitudinal relaxation 
time, as compared to the mean-field analysis. 



VI. CONCLUSIONS 

In this paper we have presented a comprehensive de- 
scription of the critical dynamics at structurally in- 
commensurate phase transitions. Our starting point 
was the time-dependent, relaxational Ginzburg-Landau 
model with 0(2) symmetry. To be more general, we dis- 
cussed the 0(n)-symmetric functional. Hence, we were 
able to study the influence of the n — 1 Goldstone modes 
accurately. We used the renormalization group theory 
in order to compute the dynamical susceptibility below 
and above the critical temperature Tj to one-loop or- 
der. Thus we could venture beyond the usual mean-field 
description. As we calculated the renormalization, fac-. 
tors in the generalized minimal subtraction schcmcj£j't£l 
we could-deal with the interesting crossover scenarios 
carefully.lij 

Our findings were used to interpret experimental data 
from NMR experiments, measuring the relaxation rate. 
The relaxation rate is connected with the calculated sus- 
ceptibility via an integral over the wavevector, at fixed 



frequency. 

Above the critical temperature T/, we showed how scal- 
ing arguments lead to an identification of the dynami- 
cal critical exponent for the relaxation rate and provide 
a qualitative understanding of its temperature depen- 
dence. Then we described the crossover from the critical 
region to a high-temperature region, where fiuctuations 
should not change the classical critical exponents. Excel- 
lent agreement for both the critical exponents resulting 
from the scaling arguments and the description of the 
crossover regions with the experimental data was found. 
This led us to the conclusion that the experimental data 
should probably not be interpreted by identifying a criti- 
cal region of supposed width of lOOif, but rather through 
a crossover between the non-classical critical exponents 
and the mean-field exponents, taking place at a temper- 
ature approximately equal to T/ -f bK. This conjecture 
yields a considerably more reasonable width of the criti- 
cal region. 

Below the critical temperature, we analyzed the dy- 
namical susceptibility calculated to one-loop order in 
the renormalization group theory in considerable detail. 
The coupling of the OP modes was considered explicitly. 
We thus gained new insight into the influence of Gold- 
stone modes on the structure of the susceptibility and 
its temperature dependence. As a result we found that 
the relaxation rate of the phason fluctuations becomes 
temperature-dependent. This temperature dependence 
disappears in the two limits when either the temperature 
approaches the critical temperature T/, or the temper- 
ature is very low. For the amplitudon modes the influ- 
ence of the Goldstone modes is more subtle. We summa- 
rized the effect in a wavevector- and frequency-dependent 
"mass" and showed that this can be interpreted as a 
bending-down of the temperature-dependent relaxation 
time as compared to a hypothetic situation where no 
Goldstone modes are present. All experimental find- 
ings are well understood treating the OP modes beyond 
their mean-field description. As reported from the anal- 
ysis ofoLltrasonic attenuation experiments for Rb2ZnCl4 
beforeJj'E3 no phason gap had to be introduced. Recently 
however, the direct observation of a "phason gm" has 
been reported for a molecular compound (BCPS)o This 
"phason gap" was observed in inelastic neutron scatter- 
ing experiments for very high frequencies. Again, the low 
frequency dynamics probed by NMR did not reveal any 
gapJia 

Thus, an interesting application of the 0(2)-symmetric 
model is presented here, in terms of a crossover descrip- 
tion and a discussion of the full k and u dependence 
of the susceptibility calculated to one-loop order. We 
found very good agreement with experimental data. Be- 
sides the precise calculation of critical exponents as one 
strength of the renormalization group theory, also de- 
tailed analysis of crossover scenarios and the effect of the 
anharmonic coupling of modes is possible. We want to 
stress how successfully the results of the renormalization 
group theory can be applied to specific experimental find- 
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ings. In addition, we emphasize that the choice of two ACKNOWLEDGMENTS 
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APPENDIX A: 



In this appendix we list analytical results for the two-point vertex functions and their singularities (1/e poles) in the 
generalized minimal subtraction scheme, following from the dynamical functional Below T/, the corresponding 

zero- and one-loop contributions are stated. For the explicit calculation of the integrals in the generalized minimal 
subtraction scheme we refer to Ref. All integrations over internal frequencies have already been performed by 
means of the residue theorem. We restrict ourselves to the three-dimensional case (e = 1). 
T > Tf. 

Here only the simplest one-loop graphs enter the diagrammatic expansion. 
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The diagrammatic expressions for the vertex functions are depicted in Figs. |6| and |T^. 
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To CTo-(k, w) =i(j + Ao (fc^ + mf) 
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u = 0.6695 



r; = 0.033 



z = 2.024 



TABLE I. Critical exponents-HH-ithe 0(2) symmetric 4>'^ 



field theory in three dimensions.! 
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TABLE IL Diflterent scaling laws for the longitudinal sus- 
ceptibility. The functions f{k) denote regular functions. 
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FIG. 1. The flow of the renormalized mass parameter t(1) 
vs. the scaling variable x = 1/t{1)'''^^^ for ten different values 
of r(l) < 0.1 [u(l) = 0.6, n = 2]. 



1.0 



0.8 



0.6 



0.4 



0.2 




10 

FIG . 2. Logarithmic derivative of the integral in Eq. 



(4.33) with respect to I, plotted vs. the scaling variable 
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FIG. 5. Propagators below Ti. 



FIG. 3. Relaxation time Ti vs. the deviation AT fromithe 
critical temperature Ti. The NMR experimental dataJ23 in- 
dicated by circles, are compared with the theoretical results, 
represented by the solid lines. A crossover to the mean-field 
regime starting at AT ~ IQK [line (a)] and at AT ~ Sif [line 
(b)] is considered here. 
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FIG. 4. Relaxation time T\ in a logarithmic plot vs. the 
deviation AT from-ithe critical temperature Ti. The NMR 
experimental datati (circles) are compared with the theoret- 
ical result, represented by the solid line. A crossover to the 
mean-field regime starting at AT ~ 5K is considered. The 
four labeled temperature regimes are discussed in the text. 
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FIG. 6. Vertices below T/. 
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FIG. 7. Counterterm below T/. 
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FIG. 8. The effective exponent We// = ^ g'fo'^i''"' of the flow 
of the renormahzed mass parameter m{l) vs. I for seven dif- 
ferent m(l) < 1 below Tj [u(l) = 10"^ ■ u*h ,n = 2, e = 1] 
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FIG. 10. The longitudinal susceptibility for the critical 
temperature T = Tj and seven other temperatures plotted 

vs. fc • y''^. The scaling regions listed in Table |l^ are valid for 
temperatures not too close to Ti. 
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FIG. 9. Flow of the coupling constant u{l) vs. / for five 
different m(1) below T/ [m(l) = 0.01, n = 2, e = 1]. The co- 
existence and the Heisenberg fixed point are marked on the u 
axis. 
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FIG. 12. The transverse susceptibility for the critical tem- 
perature T = Ti and seven other temperatures plotted vs. 

k ■ ^J^- The temperature dependence is only present in an 
interim region. 




FIG. 14. Real and imaginary part of the effective "mass" 
(self energy) / and its k dependence for Cb = 0.01, 1, and 100. 
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FIG. 13. The effective exponent k of with respect to 
k plotted vs. k - J ^ for the same temperatures as in Fig. |l^. 
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FIG. 15. Relaxation time T^ (phason dominated) and 
Tj" (amplitudon dominated^ vs. reduced temperature. The 
NMR experimental data,E3 again indicated by circles, are 
compared with the theoretical results, represented by the solid 
lines. As in the high-temperature phase a crossover to the 
mean-field regime starting at AT ~ 5K is considered. 
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FIG. 16. Transverse susceptibility and its diagrammatic representation. 
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FIG. 17. Longitudinal susceptibility and its diagrammatic representation. 
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